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We report the first detection of the Higgs-type amplitude mode using Bragg spectroscopy in 
a strongly interacting condensate of ultracold atoms in an optical lattice. By the comparison of 
our experimental data with a spatially resolved, time-dependent bosonic Gutzwiller calculation, we 
obtain good quantitative agreement. This allows for a clear identification of the amplitude mode, 
showing that it can be detected with full momentum resolution by going beyond the linear response 
regime. A systematic shift of the sound and amplitude modes' resonance frequencies due to the 
finite Bragg beam intensity is observed. 
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In recent years, remarkable progress has been made in 
the field of ultracold atoms, enabling the simulation of 
strongly interacting quantum systems beyond the scope 
of traditional solid state counterparts [THS]. In solid 
state systems, spectroscopic techniques such as ARPES 
and neutron scattering have been established as refer- 
ence methods for providing energy and momentum re- 
solved insight into the excitational structure of materi- 
als. Recently, spectroscopic techniques have also been 
applied successfully to ultracold atoms, such as RF spec- 
troscopy [1] , lattice shaking [51 [5j , as well as several ex- 
periments using Bragg spectroscopy [THTB]. Initially, the 
latter was performed on weakly interacting condensates 
[TUl H"2"] , then extended to strong interactions without a 
lattice [llj . In more recent experiments, these studies 
have also been extended to ultracold atoms in optical 
lattices [7HH1 IT3rf2"D] . which opens up the possibility of 
studying a number of models with strong correlations 
from condensed matter theory. Up to now, these Bragg 
spectroscopic experiments have, however, been focused 
on weakly interacting condensates [THS] or the Mott in- 
sulating (MI) [5] regime. In this letter, we investigate the 
excitational structure of a strongly interacting lattice su- 
perfluid (SF) and, for the first time, clearly identify the 
recently described amplitude mode [HI l2"TH2"4"] • 

At sufficiently high lattice depth, the atoms are well 
described by the Bose-Hubbard model (BHM) 

H = - J E ( b lh + h - c -) +E( e < -Mb + \ E 6 M & a , 

where b\ creates an atom at lattice site i, the tunneling 
between nearest neighboring sites is characterized by a 
tunneling matrix element J and the on-site energy shift 
an atom experiences at a given site in the presence of n 
other atoms is given by Un, with the interaction param- 
eter U. A local energy offset is accounted for by e; and 
/i denotes the chemical potential. The crucial parameter 
for realizing different regimes is the ratio U/J. In the 



Bogoliubov regime U/J <C 1, the gapless sound mode, 
corresponding to the excitation of Bogoliubov quasipar- 
ticles in a lattice has been investigated experimentally [7J- 
9 . Intermediate lattice depths allow for the realization 
of a strongly interacting SF beyond the realm of Bogoli- 
ubov theory, exhibiting a rich excitational structure: In 
addition to the gapless sound mode, the existence of the 
gapped 'amplitude' mode in the BHM (within the low- 
est band), generated by a physically similar mechanism 
as the Higgs boson in high energy physics [HI EH [22] , 
has been a topic of high interest in recent literature 
[21 [HH21] • However, linear response calculations in the 
perturbative limit have suggested that this mode can- 
not be addressed in a momentum-resolved fashion with 
Bragg spectroscopy [21] and there has been no clear ex- 
perimental signature in previous measurements [8]. To 
bridge the gap between existing idealized theory predic- 
tions and our experimental observations, we address a 
number of important experimental effects in our simu- 
lations: 1) the high probing beam intensity; 2) spatial 
inhomogeneities, such as the harmonic trapping poten- 
tial breaking the translational symmetry and leading to 
a broadening in /c-space; 3) strong interactions in the SF 
requiring a treatment beyond Bogoliubov theory; 4) the 
short probing pulse time leading to a broadened signal in 
w-space. Each of these effects can modify the resulting 
measurement and a comprehensive analysis has not been 
performed to date. 

In our experiment, 87 Rb atoms are cooled in a shal- 
low magnetic trap with ui — 2ir ■ (16, 16, 11) Hz, form- 
ing a Bose-Einstein condensate before a 3D cubic op- 
tical lattice with a spacing of a = 515nm is slowly 
ramped up to a final intensity of s recoil energies E r , 
as described in Ref. [7J. This transfers the atoms into 
a condensed state in the lowest band of the lattice, 
where the system is well described by the BHM and 
the s-wave interaction through the background scatter- 
ing length is parametrized by the interaction constant U. 
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Subsequently, two Bragg laser beams with a slight fre- 
quency detuning ujb but essentially the same wavelength 
A = 781.37nm (i.e. \wg\ <C c/A), lying in the x-y-plane 
of the optical lattice at a coincident angle 9b = 45°, are 
applied. This allows the atoms to undergo a two-photon 
process, in which the momentum kick an atom experi- 
ences is given by |p B | = (47r/A) sin(20s). Our specific 
experimental setup allows the system to be probed along 
the nodal direction. For brevity, all dispersion relations 
and results shown in this paper are along this line, con- 
necting the T = (0,0,0) and M = (1,1,0) points in the 
first Brillouin zone (BZ). 

Within a classical treatment of the laser field and 
using the correspondence principle, the effect of the 
time-dependent Bragg field on the atoms is theoreti- 
cally described by the single particle operator B(t) = 
v ^g-twsipt^ _|_ e lulBt p pB ), corresponding to a propagat- 
ing sinusoidal potential with wave vector |p B |, where V 
denotes the Bragg intensity and we use units of h = 1. 
In free space the operator p PB = S p ap+p B a p ac ^ s as 
a translation operator in momentum space and simply 
transfers atoms into higher momentum states p B , if en- 
ergetically allowed, where a p is the annihilation operator 
for a momentum state p. However, in the presence of an 
optical lattice interactions are intensified and the multi- 
band structure and periodicity of the BZ invalidate this 
intuitive picture: multiple scattering events are enhanced 
and may lead to the occupation of a broad distribution of 
momentum components [Fig.[ljc)-(d)]. Moreover, strong 
interactions require an analysis in terms of a renormal- 
ized quasi-particle picture. In this Letter we focus on 
the physics within the lowest band, requiring all relevant 
energy scales to be lower than the band gap. 

To incorporate the Bragg operator into our dynamic 
Gutzwiller (GW) calculation, it is transformed into Wan- 
nier space via the unitary transformation obtained from 
a band structure calculation as explicated in [25], Ap- 
pendix A. This leads to a lowest band representation 
Pp B = Si j PijH 6j , with the exact intra-band ma- 
trix elements pij treated beyond the on-site and near- 
est neighbor approximation (decaying exponentially with 
|i — j|), where i,j denote the site indices. Within bosonic 
Gutzwiller theory, the variational ansatz for the many- 
body state consists of a single tensor product of states 
at each site \i/>(t)) = n®i l<^i(*))ii which correctly recov- 
ers both the atomic limit U / J — > oo and time-dependent 
Gross-Pitaevskii theory within a coherent state descrip- 
tion for weak interactions and it becomes exact in high 
spatial dimensions. For a strongly interacting conden- 
sate in the vicinity of the Mott transition, it furthermore 
includes the physics of the effective theories by Hubcr 
et al. [14j [21] . For a given trap geometry and experi- 
mental parameters, the ground state is determined and 
subsequently time evolved in the presence of the Bragg 
beam. The equations of motion are determined by min- 
imizing the action (including the time-dependent Bragg 
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FIG. 1: (Color online) Comparison of the experimental visi- 
bility (a) and theoretically predicted energy absorption (b) at 
Pb| = using a Blackman-Harris pulse of 10ms in an op- 
tical lattice with s = 13 in a 3D trap. A maximum intensity 
V = 0.27-Er and the experimentally determined total particle 
number N tot = 5 ■ 10 4 ± 33% and u = 2tt ■ (26, 26, 21) Hz 
trapping frequency were used for (b), leading to a maximum 
central density n = 1.05. The lower peak is the sound mode, 
mainly broadened by the high intensity of the Bragg beam to 
lower frequencies. The upper peak at 650Hz is the gapped 
amplitude mode, broadened mainly by the trap. Figures 
(c,d,e) show the theoretically predicted trap broadened log- 
arithmic quasi-momentum distributions in the first Brillouin 
zone at the frequencies marked by the green lines in (b). 

operator) and are equivalent to the time evolution gen- 
erated by a set of effective local Hamiltonians, coupled 
non- linearly to the states at other sites (see [25], Ap- 
pendix B). 

The pulse shape investigated theoretically here, is the 
square pulse where the Bragg intensity V is constant 
over a fixed time interval t. This leads to the charac- 
teristic sine 2 response in frequency space, as is to be ex- 
pected from time-dependent perturbation theory and can 
be seen in Fig. [3] To minimize the oscillatory response for 
a restricted Bragg pulse time, a Blackman-Harris pulse 
[26] is used to obtain the central results shown in Fig. [I] 
both in experiment and theory. Since the energy ab- 
sorbed from the Bragg pulse leads to a depletion of the 
condensate after rethermalization, the former is mono- 
tonically related to the visibility [2] and it is useful to 
compare these two quantities. The lower peak at ~ 200Hz 
in the spectra is the trap- and intensity-broadened sound 
mode, whereas the higher peak at ~ 650Hz is the am- 
plitude mode, broadened mainly by the strong density 
dependence. 

While exposed to the Bragg lasers, atoms are con- 
tinuously transferred between different quasi-momentum 
states, with k = — > p B initially being the dominant 
transition at weak interactions. With increasing U / J, 
backscattering transitions are enhanced and at longer 
times higher order transitions also become relevant. This 
can also be seen from the physical momentum distribu- 
tion n(p) = (a p a p ), which is directly related to the quasi- 
momentum distribution, as shown in Fig.mc)-(e). In the 
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low intensity and long-time limit, Bragg spectroscopy di- 
rectly probes the dynamic structure factor. For fixed p B , 
the various quasiparticle energies can be determined from 
the strongest loss in the momentum component rt(k = 0), 
gain in n(k = p B ), energy absorption or reduction in the 
condensate fraction as a function of the frequency ljb, as 
shown in Fig. [lja),(b). At large s, an additional com- 
plication arises in experiments: since the condensate is 
strongly depleted, the time of flight images are very sim- 
ilar to those of a thermal cloud within the signal to noise 
ratio. Thus, the lattice depth is ramped down linearly 
over 10ms to s — 10 after exposure to the Bragg beams. 
Subsequently, the visibility, shown in Fig. [ija), is ex- 
tracted from the time of flight image of the equilibrated 
atoms and is monotonically related to the absorbed en- 
ergy. Determining these resonance positions for a range 
of different momenta p B leads to the dispersion rela- 
tions with the Bogoliubov, amplitude and higher gapped 
modes shown in Fig. [2] and compared with other theoreti- 
cal results. A probing beam at resonance with a collective 
mode frequency induces time- and position-dependent os- 
cillations of the density and the spatial order parameters 
= (^i)- ln a theoretical description, these excitations 
correspond to coherent states of the respective quasipar- 
ticle, i.e. the most classical excitation, and are graphi- 
cally illustrated in Fig. |2jc),(d): a coherent Bogoliubov 
excitation leads to a dominant spatial and temporal oscil- 
lation of the phase of tp\ (which becomes pure for fc— >0) 
and a density wave, whereas an excitation of the ampli- 
tude mode leads mainly to an oscillation of the amplitude 
of ipi and the density modulation is strongly suppressed. 
The oscillation of at constant density can thus be un- 
derstood as a local periodic transfer of particles between 
the condensate and the non-condensate. In contrast to 
the weakly interacting case, where the quasiparticle en- 
ergies of the different modes depend approximately lin- 
early on the density (black dotted lines in Fig. [2jb)), the 
dependence in the strongly interacting case is highly non- 
trivial. The strong dependence can be understood from 
the excitational particle and hole branches, which may 
cross each other in the Mott insulator: crossing the phase 
transition into the SF, the emerging condensate couples 
the particle/hole branches in the equations of motion, 
hybridizing these and leading to avoided mode crossings 
at the previous intersection points, as is shown by the 
blue squares in Fig. J2jb). For all k, U/J and densities in 
the SF, the sound (amplitude) mode remains the ener- 
getically lowest (second lowest) lying mode. Comparing 
our theoretical results with Bogoliubov theory [57] (black 
dashed line in Fig.J^a)), excellent agreement is obtained 
in the weakly interacting limit. At intermediate interac- 
tions s = 9 (U/J ~ 8.55) and density n = 1 shown in 
Fig. [3^d)), neither Bogoliubov theory (dotted blue line), 
nor the theory presented by Huber et al. [21] for strong 
interactions (dashed green lines) apply and deviate from 
our results. Here, the dispersion relation obtained by the 



amplitude mode 



• s=5 

• s=13 



a -a - 



0.4 



1 I 0.8 r , , 1.2 

|k| [ji/a] 



(b) X 




density 




lattice site 
position 



FIG. 2: (Color online) (a, b): Theoretical dispersion relations 
for the homogeneous system in the linear response limit at 
n = 1 (a) and density-dependence at |k| = it /a (b) in the 
SF for weak (s = 5, black circles) and strong (s = 13, blue 
squares) interactions. Corresponding Bogoliubov results [27] 
are shown as black dashed lines in (a). (c,d): Illustration of 
the order parameter for a coherent excitation of the phase 
(sound) mode (c) and the amplitude mode (d) in a homo- 
geneous condensate at k = (0.8/a, 0,0) and s — 13. The 
projection of all ipi's in the complex plane is shown by the 
black ellipses: for the sound mode (c), the oscillation is al- 
most exclusively in the tangential, for the amplitude mode 
(d) mainly in the radial (i.e. in the amplitude) direction. 

dynamic GW method (black circles in Fig. [3jd)) remains 
valid and continuously connects these two limiting theo- 
ries. 

An essential effect that has to be considered in a re- 
alistic modeling of Bragg spectroscopy is the finite in- 
tensity of the probing beam. This is particularly impor- 
tant for strong optical lattices, where the typical time 
scale 1/J grows exponentially with s. As the pulse time 
is restricted by decoherence, a increasing intensity V is 
required for strong lattices. The analysis of this effect 
requires a treatment beyond the linear response of the 
system (i.e. not contained in the dynamic structure fac- 
tor) , as shown in the spectra of the full time-dependent 
GW calculation in Fig. [3j Whereas the response in the 
limit of very small V shown in the insets of Fig. |3][a)-(c) 
is given by (5-shaped peaks as expected, there is a drastic 
non-trivial broadening of the different peaks for typical 
experimental intensities V ~ 0.1E r , shown in the respec- 
tive main figures. This indicates a breakdown of the non- 
interacting quasiparticle picture of the BHM due to the 
large V. Whereas the amplitude mode's signature is gen- 
erally stronger in the energy- and n(k = 0) than in the 
n(k = p B ) profile, the scaling of its spectral weight is non- 
linear in V (and thus beyond linear response in this large 
V ■ t regime), as shown in comparison of Figs. [3](a),(b) 
and the respective insets. 
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FIG. 3: (Color online) Energy absorption (a) and quasi-momentum density (b),(c) spectra for a square pulse at high intensity 
V — O.lEr (insets: weak intensity V = 0.005i5 r ) in the intermediate interaction s = 9 regime and for Bragg momentum 
Pb| = n /a. The resonance frequencies predicted from the maxima of the high intensity energy absorption spectra (plotted as 
gray squares in (d)) contain a systematic uncertainty quantified by the FWHM of the pulse after 10ms « 3. 26/ J indicated by 
the error bars and shaded region in (d). The comparison with the true quasiparticle energies (dashed white lines in (a-c), black 
circles in (d)) reveals significant discrepancies. For comparison in (d): The blue dotted line is the Bogoliubov result, the green 
dashed lines are the results from Ref. [2 1J for the amplitude and sound modes (uj(k) = tpoV2Une^ with tpo determined by GW). 



At high intensity V the spectra are not only broadened, 
but the supposed resonance frequencies of all modes (gray 
squares in Fig. |3jd)) are systematically shifted to lower 
frequencies with respect to the true quasiparticle ener- 
gies (indicated by dashed white lines in Fig. J3ja)-(c) and 
circles in d)), consistent with RPA [H]), i.e. the quasi- 
particle energies are renormalized by the interaction in- 
duced by V. The error bars and shaded areas in Fig.^d) 
indicate the FWHM of the energy absorption profile af- 
ter t — 10ms, quantifying the systematic uncertainty in 
the extracted energies. To the best of our knowledge, 
this has not been considered in the analysis of experi- 
mental data thus far. Two further effects accounted for 
in our calculation are the frequency broadening due to 
the finite pulse time, as well as the inhomogeneous trap- 
ping potential. A shallow trap and low filling n < 1.05, 
due to the strong density dependence of the mode fre- 
quencies, are crucial for an unambiguous identification 
of the amplitude mode. We stress that only by taking all 
these effects into account, the good quantitative agree- 
ment, shown in Fig. [TJb), between theory and experi- 
ment in the spectra is achieved. In |25j . we point out the 
underlying connection to lattice amplitude modulation, 
and furthermore perform a time-dependent calculation 
for the experiments [HE] in 3D, finding good agreement 
in the absorption peak frequency at different s. 

In conclusion, we have experimentally observed the 
gapped amplitude mode of the BHM in the strongly inter- 
acting superfluid regime using Bragg spectroscopy. Good 
quantitative agreement between the experimental visi- 
bility and the theoretically predicted energy absorption 
from a time-dependent bosonic Gutzwiller calculation is 
found, but only when taking the full spatial trap pro- 
file, finite pulse time and high intensity of the probing 
beam into account. This shows that Bragg spectroscopy 
is a suitable method for probing not only the quasipar- 
ticle structure of Bogoliubov mode with full momentum 
resolution, but also of the more exotic collective ampli- 
tude mode excitation. For a clear signal of the latter in 
a strongly interacting SF, a shallow trap on the experi- 



mental side and a theoretical treatment beyond the per- 
turbative linear response regime are essential. Whereas 
a finite Bragg beam intensity is vital for a clear spec- 
troscopic response of the amplitude mode, it leads to a 
renormalization of the sound- and amplitude mode reso- 
nance energies, which has to be accounted for in a quan- 
titative comparison of experiment and theory. 
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Appendix A: The Bragg Operator 



lattice sites, we define the function 



In this section we derive the Wannier representation 
of the single particle operator associated with the Bragg 
beam. This operator B(t), which effectively describes the 
virtual absorption and stimulated emission of a photon 
via an intermediate level is most naturally expressed in 
physical momentum space, where it acts as a superposi- 
tion of two translation operators p pB and p PB with time- 
dependent complex phases. For the simulation within the 
frame of the lowest band Bose-Hubbard model, we seek 
a representation of this operator, which goes beyond the 
lowest order description of a moving wave for the local 
on-site energies. Note that whereas the Bragg operator 
B(t) = \ (e~ ioJBt p PB + e iuJBt p pB ) is not separable within 
the different spatial dimensions, the two constituents p PB 
and p PB are very well separable. We therefore first focus 
on expressing p PB in the Wannier basis for the ID case 

as Ppb, = Si j Pi-j^i^j > before composing the Bragg and 
lattice amplitude modulation operators for the full 3D 
case in terms of these matrix elements. 

One starts with the single particle basis transformation 
between the Bloch and physical momentum states 



\k. a) — 



E 



ct' k) \p = 2n qi +k) 



(1) 



where qi — n/a is the lattice momentum and a the 
band index. The coefficients c^"'^ can be explicitly ob- 
tained by diagonalizing the lattice Hamiltonian in quasi- 
momentum space. Choosing the appropriate normaliza- 
tion, these coefficients constitute a unitary matrix for a 
parametrically fixed quasi-momentum k 



n 

V- (a,fc)* (a, 



(2) 
(3) 



Furthermore it is useful to define the functions N(p) = 
|^ and K(p) = p — 2qiN(p) for the transformation 
from physical momentum to Wannier space, which allows 
the transformation between Bloch and true momentum 
states to be expressed as 



b (a) f - 

°K(p) ~ 



E 



(a,K(p)) f 
C N(p+2qim) a p+2q l m- 



(4) 



As required by symmetry, this expression is manifestly 
invariant under p i— > p + 2mqi with meZ. 

For the basis transformation and subsequent projec- 
tion onto the lowest band (i.e. considering only terms 
with a = 0) of the lattice along the dimension d with 



M(|i-jU) = 



L(<9 ^ 
p 



"N(p+q) 



-N(p) 



(5) 

which fulfills the relation M (-|i - j|, q) = M* (|i - j|, q). 
The full 3D operator can then conveniently be expressed 
in terms of this function as 



Pi = ei ° PB ' i II M (0 - J) ■ e ^ Pb • e*), (6) 

d=l 



where i,j are the 3D lattice site vectors, containing 
the integer site numbering of the cubic lattice along each 
dimension and is the unit vector along the d-th di- 
mension. Note that these matrix elements do not only 
depend on the Bragg momentum p B , but also implicitly 
on the lattice depth s. 

A short calculation explicitly verifies the property 
Pi.jiVB = 0) = 5ij along one dimension, i.e. different 
layers are not coupled by the Bragg process if p B lies 
within a plane parallel to these layers. To first order, the 
Bragg operator locally corresponds to a sinusoidal poten- 
tial (shift of effective chemical potential) moving along 
the direction p B , but higher orders also give rise to near- 
est neighbor and longer range hopping. Usually these 
higher order terms are neglected, but their significance 
increases in certain regimes, such as when p B approaches 
2qi (see Appendix C). In Fig. 1 the modulus of the ma- 
trix elements pij is plotted as a function of the discrete 
lattice distance |i — j| for different lattice depths s and 
Bragg momenta pb- 




FIG. 4: Dependence of complex matrix elements on the dis- 
crete lattice distance |i— j| (subfigure a: modulus, subfigure b: 
arg). The values for q = qi/ y/2 and s = 13 are relevant for the 
Bragg experiment addressed in the main text, while the value 
of q = 2qi is relevant for lattice modulation spectroscopy, dis- 
cussed in Appendix C. The behavior of the complex angle of 
M(i— j) is determined by q and alternates in a regular fashion 
for q = 2qi (here both curves lie on top of one another) in con- 
trast to the case of q = qi/y/2. The long-distance exponential 
decay constant is primarily governed by the lattice depth s. 
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Appendix B. Time Evolution within the Dynamic 
Gutzwiller Approach 

The variational ansatz within the bosonic Gutzwiller 
approach for the state on the total lattice is the direct 
product state \ip{t)) = Yl^ \4>i{t))i with the most gen- 
eral form of a local lattice state. This can be expanded 
in any local basis, such as the local Wannier Fock ba- 
sis \4>i(t))i = J2 n fn (t)\n)i, where f^\t) are now time- 
dependent fields within the dynamic approach. As the 
action 

s = J dt [-i(i,(t)\(d t m)))+m)\n BS +B(t)\m)} 

(7) 

becomes extremal for the physical path and the en- 
ergy function is a continuous function of the amplitudes 
fn (t)t the variation with respect to any of the fields 
f n ' (t) leads to the equations of motion. The dynamics 
defined by these equations for f$ (t) is fully equivalent to 
a physically more intuitive formulation, where the time 
evolution of the local state (at site i) is generated 

by the effective local Hamiltonian 

n i = ^b\blb i b i -[fi-B i - 1 (t)}blb i ^v*(t)b i ~r ] (t)bl (8) 

Within this picture, all local on-site operators (or sums 
thereof) are accounted for exactly, whereas the non-local 
terms of any single particle operator are decoupled as 

^ EAlW^J + b l(bi) <&?><&,»- (9) 

In our specific case this leads to the time-dependent com- 
plex parameters rfr(t) = JEjen(i)^i ~ Ej^i-^ijV'j in 
Eq. ([8]), where ipi(t) = (b{) is the local order parameter, 
directly related to the Wannier space representation of 
the condensate state and n(i) denoting the set of all sites, 
which are nearest neighbors of sites of i. The parameter 
r7i (t) thus accounts for both the non-linear nearest neigh- 
bor coupling through the hopping term in % B h, as well 
as the long range coupling to all other sites in the plane, 
mediated by B(t). 

By choosing sufficiently small time steps during the 
evolution, the numerical error is well controlled. Fur- 
thermore, various constants of motion (e.g. the total 
particle number, the energy for time independent Hamil- 
tonians, the quasi-momentum for translationally invari- 
ant lattice potentials, etc.) are also conserved during the 
dynamic Gutzwiller time evolution. Numerically, realis- 
tically large (i.e. of the order 10 6 sites), strongly inter- 
acting non-equilibrium 3D lattice systems can be treated 
for all experimentally relevant time scales within this ap- 
proach. 



At every point in time physical quantities, such 
as the total energy (evaluated in the pure BHM 
Hamiltonian) E(t) — {^{t)\hL- Bn \'^{t)) , the conden- 
sate fraction or the quasi-momentum distribution 
n(k,t) — {ip(t)\b\b k \4>{t)) are directly accessible from 
the many-particle state. The latter is easily related 
to the physical momentum distribution by (apa p ) 

\ - „{a,K(p))* (a' ,K{p)) ,,{a) t („') . ■ - r 

l^ a ,a' C N(p) C N(p) \°K(p) °K(p) / &nd consls t& Ot 

scaled replicas in the higher Brillouin zones within the 
lowest band approximation. 



Appendix C. Lattice Modulation Spectroscopy 

To draw a connection to the lattice modulation exper- 
iments by Stdferle et al. [5,6], we furthermore performed 
a time-dependent Gutzwiller calculation based on their 
experimental parameters. Usually, lattice modulation is 
understood as a periodic modulation of the nearest neigh- 
bor hopping element J, whereas Bragg spectroscopy is 
thought of as a spatial sinusoidal wave in the on-site en- 
ergies, moving through the lattice at a velocity Both 
of these approximations are only true up to lowest order. 
Here we shortly derive the exact description of these two 
processes (as used in our numerical implementation) and 
also point out a fundamental connection between lattice 
amplitude modulation and Bragg spectroscopy. 

It is useful to express the optical lattice potential op- 
erator in terms of the momentum shift operator p p (up 
to a constant) as V lat — ^^{p\ qi + P-2q t )- During lattice 
modulation spectroscopy, the lattice amplitude is tempo- 
rally modulated with a frequency u m and the perturbing 
operator becomes 

V mod = \ cosK^) {p\ qi + pl 2qi ). (10) 

Decomposing the cosine into its two complex exponen- 
tial constituents, this form of the lattice amplitude mod- 
ulation operator reveals that it corresponds exactly to 
the sum of two Bragg processes with opposite frequencies 
ljb = ±w m at fixed momenta p = ±2g; in the direction of 
the modulated lattice laser beam. However, in contrast 
to Bragg spectroscopy at q ^ 2mqi with m £ Z, where 
the most relevant lowest order term is the offset of the 
local potential, this term is irrelevant for lattice ampli- 
tude modulation (at q = 2qi), where this potential offset 
is equal for every site and thus only leads to an irrelevant 
global phase. Hence, the physically relevant lowest order 
term for lattice amplitude modulation is the temporally 
modulated coupling to the nearest neighboring sites, i.e. 
a modulation of J. 
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FIG. 2. Absorbed energy in the condensate during lattice amplitude modulation for different lattice amplitudes 
s = 8, 10, 12. Note that the resonance frequency of the lowest mode (amplitude mode at k = 0) decreases with 
increasing lattice intensity (i.e. increasing interaction strength) in accordance with the experimental observation [5]. 
The insets show the ground state density profile n; (green dots) and the corresponding order parameter distribution 
ipi (black crosses) along the x-direction through the center of the 3D trap, determined for the respective experimental 
parameters. 
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The absorbed energy, calculated within a time- 
dependent bosonic Gutzwiller calculation, is shown in 
Fig. 2 as a function of the modulation frequency and 
the exposure time. The parameters for this simulation 
were chosen in accordance with the experiment [5], i.e. 
87 Rb atoms in a 826nm optical lattice were exposed to 
a modulation of 20% of the lattice depth and the maxi- 
mum shown time of 30ms corresponds to the experimen- 
tal exposure time. We considered three lattice inten- 
sities s = 8, 10, 12 corresponding to different ratios of 
U/J, all on the condensate side of the superfluid-Mott 
insulator transition. The initial density distribution was 
chosen, such that the central density agreed with the 
central density calculated for the actual experimental 
situation (insets in Fig. 2), with a total particle num- 
ber of N = 1.5 x 10 5 in a lattice with an underlying 
uj = 2ir ■ (18, 20, 22) Hz harmonic trap and translational 
invariance was assumed in the z direction. The frequen- 
cies of the absorption peaks in our simulation agree well 
with the position of the observed maxima for the differ- 
ent lattice intensities. The experimentally observed peak 
widths are, however, considerably larger than in our sim- 
ulations, an effect which may be attributed to a finite 
temperature in experiment, as opposed to T — in the 
simulation. Moreover, an increase in the harmonic trap- 
ping frequency due to the lattice beams, which was not 
specified in [5] and accounted for in our simulation, may 
lead to a higher particle density at the center and hence 
also a greater inhomogeneity in the density profile, which 
would lead to a broadening of the absorption peaks. Due 



to this discrepancy in the peak width, we cannot claim to 
obtain real quantitative agreement with the experimental 
data, as for our Bragg experiment. The good correspon- 
dence in the peak positions does however indicate, that 
the lowest peak in the 3D measurement [5] can be asso- 
ciated with the amplitude mode at k = 0. 

We found that with lattice amplitude modulation spec- 
troscopy the Bogoliubov sound mode cannot be excited, 
also when going beyond the linear response regime. Thus 
all peaks observed in the strongly interacting condensate 
can be attributed to the amplitude and higher gapped 
modes at k — and we furthermore point out that 
the larger width of the peaks cannot directly be at- 
tributed to the condensate phase. The peak ranging from 
v = 0.9kHz to 1.9kHz for different lattice intensities in 
our simulation can be identified with the amplitude mode 
at k = 0, but seems to overlap and merges with the next 
mode in the experimental data. When crossing the phase 
transition into the Mott insulator, combined particle- 
hole modes appear at an energetically similar position as 
the amplitude mode. However, the quasi-particle struc- 
ture fundamentally changes when crossing the superfluid- 
Mott insulator transition and we emphasize that the am- 
plitude mode cannot simply be seen as the [/-mode in the 
condensate. Whereas the phase response of a coherently 
excited amplitude mode state vanishes when approach- 
ing the transition from the superfluid side, this is large 
everywhere in the Mott insulator for both the particle 
and hole branches. 



